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Abstract 

We use the quasiparticle random phase approximation with a few Skyrme density functionals 
to calculate strength functions in the J'^ = 0^,1^, and 2+ channels for even Ca, Ni, and Sn 
isotopes, from the proton drip line to the neutron drip line. We show where and how low-lying 
strength begins to appear as N increases. We also exhibit partial energy-weighted sums of the 
transition strength as functions of N for all nuclei calculated, and transition densities for many of 
the interesting peaks. 

We find that low-energy strength increases with N in all multipoles, but with distinctive features 
in each. The low-lying 0"^ strength near the neutron at large N barely involves protons at all, with 
the strength coming primarily from a single two-quasineutron configuration with very large spatial 
extent. The low- lying 1~ strength is different, with protons contributing to the transition density 
in the nuclear interior together with neutrons at large radii. The low-lying 2"'" transition strength 
goes largely to more localized states. The three Skyrme interactions we test produce similar 
results, differing most significantly in their predictions for the location of the neutron drip line, 
the boundaries of deformed regions, energies of and transition strengths to the lowest 2"*" states 
between closed shells, and isovector energy-weighted sum rules. 

PACS numbers: 21.10.Pc, 21.60.Jz 



1 



I. INTRODUCTION 



Very unstable nuclei with unusual ratios of neutron number (A^) to proton number (Z) are 
currently the subject of intense theoretical and experimental work. A particularly important 
issue is the effect on nuclear properties of a large neutron excess and the accompanying low- 
lying continuum of excited states. In ground states these conditions sometimes lead to halos, 
thick skins, and changes in magic numbers. Less is known about excited states, which we 
study here by calculating transition strengths and densities. 

Most of the work on excited states is in light neutron-rich nuclei; a low-energy strength- 
function peak has been studied for years. Recent work [1-3] has focused on enhanced 1~ 
strength in nuclei like ^^Li or in neutron-rich oxygen isotopes [4-12]. But in heavier nuclei, 
calculations are more spotty. Several groups have calculated the isovector (IV) 1^ strength 
functions, or photoabsorption cross sections, in the Ca isotopes [6, 11, 13-15]. References 
[6, 15] predict a small peak at around E ~ 9 MeV in ^"Ca. The first paper also calculates 
the strength function in '^^Ca, finding no low-energy peak, but Ref. [14], a calculation based 
on density-functional theory, finds a very broad bump ai E = 6-10 MeV. It is not clear, 
therefore, at what value of N the low-energy strength is first discernible. 

The method we use here, the quasiparticle random phase approximation (QRPA), has 
been applied mainly to selected 1~ excitations in particular isotopes of O. Ca, and Ni; one 
example is Ref. [11], which asserts that the low-energy excitations rcficct strong particle- 
particle ( "dineutron" ) correlations. References [16, 17] show extensive and systematic QRPA 
calculations of photoabsorption cross sections, but focus primarily on radiative neutron- 
capture cross sections. The relativistic QRPA [9, 18-20] has been applied more with an eye 
to nuclear structure issues (for a recent review, see Ref. [21]), asserting, for example, that 
low-lying isoscalar (IS) 1~ states exhibit toroidal fiow. References [22, 23] go beyond RPA, 
altering details of predicted low-lying strength. 

Other channels have been discussed less than the IV 1~. The authors of Ref. [24] calculate 
O"*" strength functions in Ca isotopes up to the neutron drip line, finding an enhancement 
at low energy in ^°Ca. Reference [25] examines transition strengths in the 0+, 1^, and 2+ 
channels near the neutron drip line, finding, for example, a sharp 2+ peak near ii^ = in 
ii^Ni. 

The above represents a considerable amount of work on excitations in exotic nuclei. Its 
range and results can be summarized, roughly, as follows: 1) Most of the calculations are for 
spherical nuclei, and the IV 1^ channel has been studied the most. 2) A low-energy peak 
grows as N increases. 3) The low-energy excitations do not seem collective in the 1~ channel 
(e.g. Ref. [26]), i.e. they do not contain more than a few particle-hole or two-quasiparticle 
configurations. 4) The corresponding transition densities have long neutron tails extending 
beyond the bulk nuclear radius. 5) The neutron skin thickness is correlated with the amount 
of low-energy strength [27]. 6) Calculations reproduce measurements of the amount of low- 
energy strength fairly well if they include many-particle many-hole correlations [22, 23]. 

Most of this work has focused on nuclei right at the drip line and, as already noted, has 
emphasized the IV 1~ channel. It is still not known whether theory makes similar predictions 
for other neutron-rich nuclei, or to what extent the properties of the low-energy peaks vary 
from one multipole to another. The purpose of this paper, therefore, is to calculate strength 
functions in as many nuclei as possible, in all multipoles up to 2+, and investigate how the 
properties of the strength functions vary as N and the multipolarity change. We use the 
QRPA with Skyrme density functionals (which we refer to as interactions from here on) and 



2 



volume-type delta pairing interactions. 

Currently, the QRPA is the most sophisticated method that can be used systematically 
for resonances of all energies in heavy nuclei without assuming that part of the nucleus forms 
an inert core. Here we will display results of QRPA calculations for Ca, Ni, and Sn isotopes 
from the proton drip line to the neutron drip line, describing the corresponding strength 
functions and transition densities in detail. We will also examine whether results depend on 
which Skyrme interaction we use, so that we can identify important measurements/nuclei 
that will help distinguish among them. 

Section II below describes our method. Section III shows strength functions in the three 
isotopic chains and looks for {N, Z)-dependent trends in their behavior. Section IV discusses 
transition densities to states in low-lying peaks. Section V is a conclusion. 

II. CALCULATION 

Our approach to the QRPA is described in detail in Ref. [28], and we discuss it here only 
briefly. The starting point is a coordinate-space Hartree-Fock-Bogoliubov (HFB) calculation, 

carried out with the code used in Refs. [29, 30]. The code outputs quasiparticle wave 
functions, with spherical symmetry assumed. It takes the continuum into account by placing 
the nucleus in a 20-fm box. We use as large a single-particle space as possible and solve the 
HFB equations accurately, so as to ease the removal of spurious states from the QRPA and 
satisfy sum rules. 

We diagonalize the HFB density matrix to obtain a set of canonical states — as many 
as there are quasiparticle states — that make up our QRPA basis. The QRPA-Hamiltonian 
matrix then contains off-diagonal elements from the quasiparticle one-body Hamiltonian 
as well as the matrix elements of the residual two-body interaction. The two-body matrix 
elements are functional derivatives of the energy with respect to the density and pairing ten- 
sor. After diagonalizing the QRPA Hamiltonian, we use the eigenvectors to obtain strength 
functions [31] and transition densities (see Appendix). 

We cannot handle an infinite number of canonical quasiparticle states, and so truncate 
our basis in the following way: For N < 82 the maximum single-particle angular momentum 
imax is 15/2, and for > 84 jmax is 21/2. If the ground state has a finite pairing gap in 
the HFB calculation we truncate the canonical quasiparticle basis in the QRPA by omitting 
states with very small occupation probability vf. In the J'" — 0"^ and 1~ channel we keep 
states with > lO"-*^^, and in the J'^ = 2+ channel we keep those with vf > 10~^. If the 
mean field has no pairing we use the Hartree-Fock energy to truncate; in the J'^ = O"*" and 
1~ channels we keep single-particle states with energy Si < 150 MeV, and in the J'^ = 2+ 
channel we keep those with < 100 MeV. In the HFB calculations we always cut off the 
quasiparticle spectrum at 200 MeV. 

The largest QRPA-Hamiltonian matrix that we encounter here (in the 1~ channel for 
Sn near the neutron drip line) is about 14,000, and the minimum is about 2,000 (in the 
O"*" channel for Ca near the proton drip line).^ The separation of spurious excitations from 
physical states is as complete as in Ref. [28] in most of cases. 

We use 3 Skyrme interactions: SkM*, SLy4, and SkP. Our pairing interaction has zero 



^ In Ref. [28], cut-oflFs of vf = 10 and = 150 MeV were used for the 2+ channel, and then the maximum 
dimension of the Hamiltonian matrix was about 20,000. 
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range with no density dependence (this is "vohime pairing"). For reasons we do not under- 
stand, we cannot obtain solutions of the HFB equation with SkP for Ca and Ni isotopes. 
The other interactions give solutions that are unstable against QRPA quadrupole vibrations 
for isotopes with: = 72 in Ni (SkM*); = 104 - 112 (SkM*) and N ^ 60 - 64, 96- 110 
(SLy4) in Sn, and so we do not discuss these nuclei. 



III. STRENGTH FUNCTIONS 

Strength functions are important because they capture information about all excited 
states and can often be measured. We define the strength functions through the transition 
operators 
J"" = 0+: 

i=l 

AT ^ 7 ^ 

i=l i=l 



J"" = 1- 



i=l 

N ^ 7 ^ 

^Im = '-^Y.'-^y,m{^,) -''—Y^nY.Mm, (2) 

i=i i=i 

J"" = 2+: 

i=l 

N ^ 7 ^ 

Fll = ^ E ^-mW -'-jY.'" ^^M(n.), (3) 

i=l i=l 

where A is total-nucleon number, and i is a single-particle index. The upper limit for Z (N) 
on a summation sign means that the sum is over protons (neutrons). 

In the next few sections we show the strength functions, with spurious states removed, 
produced by SkM*, and discuss those produced by the other Skyrme interactions when ap- 
propriate. We then display energy-weighted sums (up to several different energies) produced 
by all our interactions, in all isotopes. 



A. The 0+ channels 

Figures 1-3 show the IS and IV strength functions in the J'^ — 0+ channels. In Ca 
(Fig. 1), most of the functions are smooth except for small low-energy spikes in several 
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nuclei and the large low-energy peaks in ' Ca. The IS giant resonances peak at E = 20 
MeV for ^^^^^Ca and gradually move to 15 MeV by ^^Ca. As approaches the neutron 
drip line, a peak near E = grows, becoming quite large at the drip line. The IV giant 
resonance is broad in the light isotopes and develops a low-energy component as N increases, 
apparently increasing the summed strength. The larger peak varies only slightly in energy 
over the same range of N. Near the drip line the low-energy IV strength develops a peak 
that exactly mirrors the IS peak, indicating that the strength in both channels is produced 
solely by neutrons. This distinct low-energy bump seems to appear first at 10 MeV in ^°Ca. 
This nucleus is not as short-lived as the isotopes near the drip line; perhaps the appearance 
of the low-energy mode is experimentally testable. The same nucleus marks the beginning 
of the low energy peak when we use SLy4. 

The N dependence of the strength function of Ni (Fig. 2) is similar in many ways to that 
of Ca; the energy of the IS giant-resonance energy, for example, falls in a similar way past 
N — 40. A large low-energy peak again develops near the neutron-drip line, though not 
as dramatically as in Ca, and with the complication that the strength shrinks after ^'^Ni. 
This peak first is noticeable in both the IS and IV channel at 6 MeV in ®°Ni, here and 
in calculations with SLy4. This example and that in Ca suggest more generally that the 
nucleus with 2 neutrons outside a closed shell will mark the appearance of the low-energy 
peak. 

In Sn (Fig. 3), the peak near zero energy around the drip line is much smaller than in 
Ca and Ni. The bump first appears in ^^^Sn a.t E = 8 MeV (violating by two nucleons the 
conjecture made a few lines up) in both the IS and IV channels. This threshold nucleus is 
the same with other Skyrme interactions. 

Although we can make these kinds of comparisons in words, we would like to display the 
results of both (or all three where possible) Skyrme interactions in a more quantitative and 
graphic way. To limit the number of our figures, we do so by showing predictions for the sum 
of energy weighted strength up to 4 successively higher energies. We define the "partial" 
sum, up to an energy E as 

We^ dE'E'S{E'), (4) 

where S{E') is a strength function^ (see Eq. (1) of Ref. [28]). Figure 4 shows We with 
E = 10, 20, 30, and 100 MeV for the IS 0^ channel in Ca. All Skyrme interactions show 
similar A"-dependence and £^-dependence; the difference between the £^ = 30 and 100 MeV 
curves is comparable with or smaller than the differences among the other curves (separated 
by 10 MeV), reflecting the eventual vanishing of the strength as E increases. The E = 20 
MeV curve increases more rapidly with A^ than the others in the region A^ > 30, indicating 
that the strength around £^ = 20 MeV is shifted down at large A". 

The differences among the interactions are most apparent in the location of the neutron 
drip line, a static property. The differences in strength functions, at least in this channel, 
are for the most part minor (though large enough to make significant differences in, e.g, the 
compression modulus). 



^ If the strength function has a large peak near E = 0, Eq. (4) overestimates the energy- weighted sum, 
because the tail in the negative energy region is not included. This problem occurs near the neutron drip 
line, and the energy-weighted sum is then overestimated by 10 % at most. We did not calculate We from 
negative E', because then WioMeV is sometimes negative. 
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Figure 5 shows the partial sums in the IV channeL Here the two interactions do yield 
somewhat different results. Though the A^-dependence is similar, Sk;M* produces more 
integrated strength than SLy4. To see how the total energy-weighted strength can differ 
from one interaction to another, we examine the IV energy-weighted sum rule (EWSR) 

5^ ^£;,|(fc|F}](,|0)r = 5ew(J) + CMJl (5) 

k M 

Cc^^ih+t2) + ^{hXi+t2X2), (8) 

^c(^) = E / \'^fjMir)\'p-oir)pl{r), (9) 



/oo(r-) = r', (10) 
fjM{r) = r'YjMm, J>1. (11) 

Here Ek is the energy of the excited state |A;), |0) is the ground state, m is the nucleon mass, 
and {r^)p {{r'^)n) is the mean square proton (neutron) ground-state radius. The sum over 
states k does not include spurious states. The second term CcIc{J) in Eq. (5), with p^{r) 
and p^ir) of Eq. (9) the proton and neutron ground-state densities, is a correction arising 
from the momentum-dependent terms of the Skyrme interaction. The Skyrme parameters 
to,ti, ■ ■ ■ are defined in Refs. [32] (SkM*), [33] (SLy4), and [29] (SkP). Another term in the 
sum rule depending on the proton-neutron mixed density is always zero in our work. 

In Tab. I we show 5'ew(0), Cc, /c(0), and EWSR (the full energy- weighted sum) for ®°Ca. 
The interactions reproduce ground-state observables about equally well, so the difference in 
the EWSR, as the tabic shows, comes from the term involving Cc- The degree by which 
the EWSR exceeds S-ew{J) is in line with expectations (see, e.g.. Chap. 10 of [34] and 
Refs. [33, 35]). 

The properties of Ni are similar to those of Ca. W£;=20MeV in the IS channel of Ni increases 
more rapidly than the other curves (Fig. 6). SkM* again has the larger IV sum-rule value. 
(Fig. 7). We=20McV in the IS channel of Sn is much higher than We=wmcV (Fig. 8), reflecting 
the lower energy of the giant resonance. Figure 9 shows that the IV sum rules follow the 
same order as in the other isotopes, with SkP shghtly larger than SLy4. 
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TABLE I: Components of the EWSR of IV 0+ channel of 6°Ca. 



parameter set 


'S'ew(O) 
(e^MeVfm^) 


(MeVfrn^) 


W) 

(eVfm) 


EWSR 
(e^MeVfm^) 


SkM* 


15587 


68.74 


58.01 


19575 


SLy4 


15722 


32.47 


57.90 


17602 


SkP 




44.64 
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FIG. 1: IS (sohd) and IV (dashed) 0"*" strength functions for even Ca isotopes (SkM*). 
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FIG. 2: The same as Fig. 1 but for Ni isotopes. 
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FIG. 3: The same as Fig. 1 but for Sn isotopes. 
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FIG. 4: AT-dependence of the partial energy-weighted sums in the IS 0"*" channel of Ca. The upper 
limits of the sums are 10 MeV (plus symbol), 20 MeV (x), 30 MeV (asterisk), and 100 MeV (open 
square). 
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FIG. 5: The same as Fig. 4 but for the IV channel. The open squares now correspond to an upper 
limit of 150 MeV. 
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FIG. 6: The same as Fig. 4 but for Ni. 
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FIG. 7: The same as Fig. 5 but for Ni. 
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FIG. 8: The same as Fig. 4 but for Sn. 
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B. The 1 channels 



Having analyzed the 0+ strength functions in detail, we will, for the most part, let the 
higher-multipole figures speak for themselves. Just a few remarks about the 1~ channel: The 
IS strength functions in Ca appear in Fig. 10. Here, unlike in the 0+ channel, there is a clear 
low-energy peak in all nuclei, increasing in size^ after ^^Ca. (The IS strength functions of 
72,74,76Qg^ are multiplied by 0.1 in the figure.) In ^^Ca 45 % of the energy-weighted strength 
is below 10 MeV. The IV strength functions also have increasing low-energy strength at 
large N. The small low-energy bump in ^'^Ca is produced by SLy4 as well as SkM*. 

We see similar iV-dependence in Ni (see Fig. 11), though after ^°Ni the low-energy IS 
peak shrinks. The threshold for rapid growth with both Skyrme interactions is the closed- 
shell nucleus ^^Ni. In Sn, once again, the low-energy IS peak grows with N, though more 
smoothly than in the other isotopes (see Fig. 12). In particular, there is apparently no 
sudden increase at the closed-shell-plus-two nucleus ^^^Sn. Recently, Ref. [36] reported a 
clear IV bump around 10 MeV in that nucleus. The bumps in our calculated spectra do not 
seem to really set in until higher N. 

The 1~ We tell a shghtly different story (Figs. 13-18). The partial sums in the IS channel 
of Sn (Fig. 17) show a noticeable kink &t N — 82, due clearly to an increase in low-energy 
strength. As remarked above, the increase is not obvious in the strength-function figures. 
The choice of interaction makes very little difference in the We except in the blank regions 
(indicating quadrupole deformation) and near the neutron drip line. 

According to the analytical argument in the previous subsection, if the IV EWSR depends 
on the interaction in one multipole, then it should in all multipoles. As Figs. 14, 16, and 18 
show, this is indeed the case. 

We have not commented on the question of whether protons are involved in the transitions 
to low-lying states in this channel because the IS and IV operators do not have the same 
radial dependence. We will address the issue in later when we display transition densities. 



^ The low-energy peaks of "'^Ca are so sharp because we smooth the strength functions with a smaUer 
width at bound states than in the continuum (see Ref. [28]). 
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FIG. 10: IS (solid, scale on left) and IV (dashed, scale on right) 1 strength functions for even Ca 
isotopes (SkM*). The IS strength functions in ^^'^^'^^Ca are reduced by a factor of 10. 



14 




FIG. 11: IS (solid, scale on left) and IV (dashed, scale on right) 1 strength functions for even Ni 
isotopes (SkM*). 
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FIG. 12: The same as Fig. 11 but for Sn isotopes. 
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FIG. 13: The same as Fig. 4 but for the 1 channel. 
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FIG. 14: The same as Fig. 5 but for the 1 channel. 
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FIG. 15: The same as Fig. 6 but for the 1 channel. 
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FIG. 16: The same as Fig. 7 but for the 1 channel. 
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FIG. 17: The same as Fig. 8 but for the 1 channel. 
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FIG. 18: The same as Fig. 9 but for 1 channel. 
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C. The 2+ channels 



The strength functions in the 2"*" channels are quahtatively different from those in the 
lower-multipole channels. In Ca, for example (Fig. 19), all isotopes^ except N — AO have 
considerable low-energy strength, in many instances because of sharp quadrupole vibrations. 
The low-energy strength grows with A^, but unlike in the O"*" channel, both neutrons and 
protons are involved in the transitions to these states, even near the drip line, because IS 
and IV distributions do not coincide. This is so for all Skyrme interactions. In fact the IS 
strength at low energies is much larger than the IV strength in all isotopes except ^^Ca. 
Even there, the IS strength to the low-energy peak with SkM* is 256 e^fm^ and the IV 
strength only 136 e^fm^. 

The 2"*" strength functions in Ni (Fig. 20) are similar, though with two differences: all the 
nickel isotopes have low-energy peaks, and for ^^^^^^Ni the IV peak is at a different energy 
than the IS peak. The lowest-lying peak seems to be nearly pure IS, like the familiar surface- 
quadrupole vibrations in stable nuclei. In the next section we will discuss the extent to which 
the states near the neutron-drip line display the same vibrational features we observe near 
stability. 

SLy4 puts the drip line at ^^Ni, much closer to stabihty, and in isotopes close to ^^Ni 
puts the IS and IV peaks at about the same location. In the lighter {A > 80) Ni isotopes, 
the lowest peak always has a little bit of IV strength mixed with the IS strength, as with 
SkM*. In s°Ni (SkM*), the IV strength in the peak at 400 keV is 5 % of the IS strength. 

Things are quite similar in Sn (Fig. 21); for ^^^"^^Sn the lowest-energy peak has both 
IS and the IV components but in ^^^~^^^Sn it has no IV component. The other interactions 
give the same results near the neutron drip line, indicating again that proton excitations are 
important in these modes everywhere. 

The predictions for We in Ca (Fig. 22) show that SLy4 evinces slightly more of a closed 
shell (a small kink in the 10-MeV curve) at = 40 than does SkM*, while showing less of a 
shell effect at = 28 (no small kink in the 20-McV curve). As was the case in the 0"^ and 
1~ channels, the IV We (Fig. 23) depends on the interaction. 

Similar differences in shell-strength among the Skyrme interactions are apparent in the 
Ni and Sn We curves for the IS channel (Figs. 24 and 26), while the difference is not obvious 
in the IV channels (Figs. 25 and 27). The IS We for Sn (Fig. 26) indicate that SkP has 
more low-energy strength (A^ > 82) than the other parameter sets. The IV curves in Sn 
show interaction differences similar to those in the 1~ channel. 



^ Reference [15] calculates the 2+ strength functions of ^°Ca without low-energy peaks; we do not know the 
reasons for the difference between that calculation and ours. 
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FIG. 19: IS (solid) and IV (dashed) 2+ strength functions in even Ca isotopes (SkM*). 
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FIG. 20: The same as Fig. 19 but for Ni isotopes. 
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FIG. 21: The same as Fig. 19 but for Sn isotopes. 
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FIG. 22: The same as Fig. 4 but for the 2+ channel. 
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FIG. 23: The same as Fig. 5 but for the 2+ channel. 
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FIG. 24: The same as Fig. 6 but for the 2+ channel. 
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FIG. 25: The same as Fig. 7 but for the 2+ channel. 
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FIG. 26: The same as Fig. 8 but for the 2+ channel. 
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FIG. 27: The same as Fig. 9 but for 2+ channel. 



25 



IV. TRANSITION DENSITIES 



Giant resonances lie in the continuum and their A^-dependence is therefore smooth. The 
low-energy peaks, on the other hand, can vary dramatically with N, especially near the drip 
line. In this section, to learn more about their collectivity and N dependence, we examine 
transition densities to representative low-lying states. 

The transition density p1^{r;k) is defined in the Appendix. Here we focus on the radial 
transition density 

pl{r;k) = r'J dnYj^MMpUr; k). (12) 
The transition amplitude depends on this quantity through the relation 

{k\FjM\0) = JdrY^ f\r)pl{r;k), (13) 

q=p,n 

where /^(r) is the radial part of the multipole operator Fjm- The radial transition density 
therefore specifies the radius at which the strength is large or small. The transition density, 
as Eq. (A. 2) shows, is proportional to the product of particle and hole wave functions and so 
is localized as long as all hole states are localized, i.e. as long as the ground state is bound. 
In this section, as before, we will be using SkM* unless we state otherwise. 



A. The 0+ channels 

Figure 28 shows the transition densities to 4 states below the giant resonances in ^^Ca with 
large strengths (the energies label the results of the QRPA calculation before smoothing). 
Proton excitations are important, particularly because of the factor in the 0+ transition 
operator. The rms radii of the protons and neutrons in the ground state are 3.40 and 3.19 
fm. The important part for the strength comes from outside the rms radius. 

This is not true of the giant resonances. Figure 29 shows transition densities in the peaks 
of the giant resonances of ^^Ca. The upper two panels represent the IS giant resonance, 
and as is expected the protons and neutrons have similar transition densities. In the lower 
panels, corresponding to the IV giant resonance, the two densities are out of phase. The 
peaks of the proton transition densities are closer to the center of the nucleus than those of 
the lower-energy states, particularly in the IS channel. 

We saw that a small bump appeared after magic numbers in neutron-rich nuclei. In 
Ca, the relevant isotope is ^°Ca. (Fig. 30) shows that the excitation is created entirely by 
neutrons. But the state is not at all collective. The largest two-quasiparticle component in 
the QRPA wave function is i^3p3/2{£ = 9.175 MeV, vl = 6.0 x 10"^) ® iy2p3/2{S = 0.880 
MeV, = 0.488), where S is the quasiparticle energy in the HFB quasiparticle basis and 
v"^ is the diagonal ground-state occupation probability (norm of the lower component of the 
quasiparticle wave function). The radial quantum numbers are assigned in order of S {—S 
for those having v'^ > 0.5). X^j^yiY^y) is the forward (backward) amplitude of the QRPA 
solution. 

The squared amplitude X^^^ — Y'^^ of the z^3p3/2 ® v^i'Viji component is 0.980. The com- 
ponent corresponds to the excitation of a neutron from a J93/2 orbit near the Fermi surface 
to another P3/2 orbit about 9 MeV above. The first orbit is strongly affected by pairing cor- 
relations. The two-quasiparticle energy is 10.055 MeV, and the energy of the excited state 
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itself is 10.005 MeV; the excitation is clearly almost unaffected by the residual interaction. 
Yet this nearly pure two-quasiparticle state, because the neutron transition density has such 
a long tail, carries 5 % of the total IS transition strength. Furthermore, protons contribute 
essentially nothing to the transition; the maximum of the X'^^ — Y^^j^ for protons is 0.3 x 10~^. 
This remarkable situation — dominance by a single two-neutron-quasiparticle configuration 
with very large strength, occurs in the O"*" channel repeatedly. 

Take, for example, the drip-line nucleus ^^Ca. Two low-lying states, ai E = 2.391 and 
2.995 MeV, together account for 25.8 % of the total IS strength and 10.2 % of the EWSR. 
Transition densities, shown in Fig. 31, are again dominated by neutrons, with outer peaks 
at 9-11 fm. The rms ground-state radii are 3.738 fm for protons and 4.615 fm for neutrons, 
so the highest (inner) peaks of ptr{i^', k) are close to the neutron rms radius. The main 
component of the excited state with E = 2.391 MeV, with a squared amplitude of 0.807, 
is i/3si/2(^ = 0.355 MeV, = 0.645) ® u4:Sy2iS = 2.092 MeV, vl = 0.012). The main 
component of the state at E = 2.995 MeV, with a squared amplitude of 0.930, is v2d^i2{£ — 
0.617 MeV, vl = 0.78) ® uSd^/iiS = 2.457 MeV, vl = 6.1 x lO"''^). Like before, these are 
nearly pure two-neutron-quasiparticle states with long tails in their product. 

Turning now to Ni near the closed N — 50 shell, we see that ^°Ni has a small bump around 
5 MeV, while ^^Ni does not (Fig. 2). The transition density in Fig. 32 again shows neutron 
dominance. This excited state has 2.8 % of the total IS strength, and its main component 
is u2d5/2{S = 0.738 MeV, vl = 0.31) u3d5/2{S = 5.088 MeV, vl = 7.1 x 10^^), with 
Xl^-Y^^ = 0.761. At the drip-line nucleus ^^Ni, the main component of the excited state at 
E — 3.452 MeV, which is one of the states forming a low-energy bump, is {i'3d3/2{£ — 1.799 
MeV, vl = 0.359)}2; it has X^^ - Y^^ = 0.818. 

In Sn the nucleus just past the closed neutron shell is ^^^Sn. Here, as we noted earlier, 
a new bump is not visible (Fig. 3), until ^^^Sn. We show the transition density to the 
corresponding state in Fig. 33. The excitation's main component is v2fij2{£' = 0.950 MeV, 
vl = 0.44) ® v?>hi2{E = 6.302 MeV, vl = 6.2 x 10"^), with Xl^ - Y^^ = 0.902. £ = 0.950 
MeV is the smallest neutron quasiparticle energy of the HFB calculation. Figure 34 shows 
the transition densities to three states in the low-energy bump in ^^^Sn. The behavior is 
similar to that in the drip-line Ca and Ni nuclei. 

We turn finally to the dependence of the transition densities on the Skyrme interaction. 
Figure 35 shows, unsurprisingly by now, that SkM* and SLy4 give similar transition densities 
to states below the giant resonances in ^^Ca. We show transition densities for some giant- 
resonance states in Fig. 36. The IS/IV character of those states is clear. 

Next we compare predictions for ^°Ca. Figure 37 shows the SLy4 transition densities 
to the states in the small bump in the low-energy continuum. The bump emerges there 
with both interactions. ^^Ca is the neutron-drip-line nucleus obtained with SLy4; Figure 38 
displays the strength functions in that nucleus. Apparently they are similar not to those 
of the SkM* drip-hne nucleus ^^Ca but rather to those of the SKM* ^^Ca (see Fig. 1). 
Clearly the quantum numbers of the quasiparticle orbits involved in the transition affect the 
transition density much more than the chemical potential. 

Figures 39 (SLy4) and 40 (SkM*) show the transition densities to the excited states in 
the shoulder at about 9 MeV in ^^Ca. The positions of the outer neutron peak are nearly 
the same in both, but the neutron transition densities for r < 4 fm have some differences. 
The proton transition density is small except in the right panel of Fig. 40. The tails of 
the transition densities are quite different from those in Fig. 31, for the nucleus with 10 
more neutrons. The two-quasiparticle components with the largest amplitude involve either 
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FIG. 28: Proton and neutron transition densities to the four states in ^^Ca in the peaks of the 
strength function {E < 17 MeV) in Fig. 1. 




FIG. 29: The same as Fig. 28 but for the giant resonances. The states in the upper panels 
correspond to the peaks of the IS giant resonance, and in the lower to the IV giant resonance. 



neutron 5^9/2 's or p3/2's, depending on the state and interaction, and the largest amplitude 
itself ranges from 0.5 to 0.99. 
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FIG. 30: Transition densities to the excited state in the small low-energy bump in the 0"'" strength 
function of ^°Ca. 
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FIG. 32: Transition densities to the excited state in the small low-energy bump in the 0"*" strength 
function of ^°Ni. 
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FIG. 33: Transition densities to the excited state in a small bump in the 0"*" strength function of 
i36Sn. 
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FIG. 34: Transition densities to the excited state in a small bump in the 0"*" strength function of 
i^esn. 
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FIG. 35: Transition densities to the 0+ excited states below the giant resonances in ^^Ca with 
relatively large strength (SLy4). 
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FIG. 39: Transition densities to the excited states in the low-energy shoulder (IS) or peak (IV) 
around E = 9 MeV of the strength function of Fig. 38 (SLy4). 
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FIG. 40: Transition densities to the excited states in the low-energy shoulders around E = 9 MeV 
and 11 MeV in the 0+ strength function of ^^Ga (SkM*). 
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B. The 1 channels 

As we have just seen, the large neutron tails in the transition density are responsible for 
enhanced low-lying strength near the drip line. The robust enhancement of low-lying IS 1" 
strength is due in large part to the factor r^, which emphasizes large radii, in the transition 
operator. 

One stark difference between these transition densities and those of the 0+ channel is the 
role played by protons, which is non-negligible here. Figure 41 shows the transition densities 
of two states in the low- lying IS peak in ^'^Ca. Protons and neutrons contribute coherently. 
The component of the state at E = 8.228 McV with the largest amplitude is v2p^i2{£ = 0.880 
MeV, vl = 0.488) ®z/lrf3/2(^ = 8.219 McV, = 0.981). The squared amplitude X^^-F^i is 
only 0.36. Two other configurations, f^3si/2®f^2p3/2 and 7r2p3/2®7rl(i3/2 have X^^— = 0.20 
and 0.12, not much smaller. And since the two-quasiparticle energy of the largest component 
is 9.099 MeV, the excited state gains more correlation energy than those in the O"*" channel. 

The proton contributions remain even at the neutron drip line. Figure 42 shows the 
transition densities to 3 excited states in the large low-energy IS peak in ^^Ca (Fig. 10). 
These states have IS strengths 9227 e^fm^ {E = 1.931 MeV), 3219 e^fm^ {E = 2.199 MeV), 
and 4856 e^fm^ {E = 3.501 MeV), representing 42 %, 15 %, and 22 % of the total strength. 
While protons play no role for r > 6 fm, they are active around r = 5 fm. The main difference 
between these states and those of ^'^Ca is that here the neutrons both extend much further 
out and contribute more inside the nucleus. It is the tail that is mainly responsible for the 
increased strength. 

Interestingly, the states in ^^Ca are mainly of two-quasiparticle character even though 
there is a coherent proton contribution. The component in the state at E = 1.931 MeV with 
the largest amplitude (X^^-F^^^ = 0.914) is u3si/2{S = 0.355 MeV, = 0.645) ®z/3p3/2(£^ = 
1.582 MeV, = 0.592 x lO"-'^). In the state at E = 2.199 MeV the component with the 
largest amplitude (X^^ - Y^^ = 0.975) is i/3p3/2 i^24/2(£ = 0.618 MeV, = 0.776), 
and at = 3.501 MeV it is 1/2^/2 1^2/^/2(8 = 2.860 MeV, vl = 0.1 x 10"^), with 
X^y — F^j^ = 0.936. The proton two-quasiparticle configuration with the largest amplitude 
has X2, - Y^^ = (0.2 - 0.3) x lO'^. 

In Ni, a low-energy peak rises suddenly at ^^''Ni (Fig. 11). The transition densities of the 
states in the peak, shown in Fig. 43, are similar to those of the neutron-rich Ca isotopes. The 
density in ^^Ni, at the drip line, is similar though larger. In Sn there is no sudden increase 
in the low-energy peak around X = 82 (see Fig. 12). Fig. 44 shows the transition densities 
to low-energy states in three isotopes around X = 82. Here the proton contributions stretch 
to larger radii than in Ca or Ni, even though they are dwarfed by the neutron tails there. 

The other Skyrme interactions give very similar results, both in the IS and IV channels 
and for all three isotope chains. 
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FIG. 41: Transition densities to the excited states in the low-energy peak of the IS 1 strength 
function of ^°Ca. 
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FIG. 42: Transition densities to the excited states in the low-energy peak of the IS 1 strength 
function of ^^Ca. 
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FIG. 43: Transition densities to the excited states in the low-energy peak of the IS 1 strength 
function of ^°Ni. 
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FIG. 44: Transition densities to the excited states in the low-energy peak of the IS 1 strength 
function of ^^'^Sn (upper panels), ^^^Sn (middle), and ^^^Sn (lower). 
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C. The 2+ channels 



Figures 19-21 showed that at the neutron drip hne the IS and IV strength functions in 
Ni and Sn have distinct low-energy peaks (though the IS strength is always larger, even at 
the low energy IV peak) while in Ca the peaks coincided. This is true, except in Ni with 
SLy4^, no matter which Skyrmc interaction we use. With SkM*, the IS strength of the 
large low-energy peak at E = 0.096 MeV in ^®Ni is 92 % of the total strength (8.9 % of 
the EWSR), and that of the peak at E ^ 0.713 MeV in ^^^Sn is 23 % of the total strength 
(2.2 % of the EWSR). In Fig. 45 we compare the transition densities of the states in these 
low-energy peaks with those in ^^Ca. As Fig. 19 suggests, though the protons play a role, 
the neutrons are proportionately more important in Ca and do not seem to correlate with 
the protons. In ^^Ni and ^''^Sn, by contrast, the neutrons and protons are completely in 
phase (as they must be if the IV strength is negligible). 

The other Skyrme interactions bear out these conclusions, though in the SLy4 calculation 
of the drip-line Ca isotope the ratio of proton transition density to neutron transition density 
is larger than that with SkM* (shown in the left panel of Fig. 45). In all nuclei near 
the neutron drip line, the states in low-energy peaks are mostly above neutron-emission 
threshold, yet the neutron tails cut off at much smaller radii than to those in the 0"^ and 1~ 
channels. In addition, the transition densities in those other channels generally have a node 
at about r = 5 fm, while the 2+ transition densities show no such feature. 

We put off for just a moment the obvious question of whether these low-lying IS peaks 
are familiar quadrupole vibrations and turn briefly to the IV channel. Figure 46 shows the 
transition densities to states in the low-energy IV strength-function peak in ^'''^Sn and ^^Ni. 
The states in ^''^Sn display remarkably IS character, despite being at the peak of the IV 
distribution. The fact that the IV distribution peaks there is a little misleading because the 
IS strength is still considerably larger, even though it is larger still at lower energies. We get 
some IV strength because the ratio of the transition density of protons to that of neutrons 
is not exactly equal to Z/N (see Eq. (3)). 

Returning to the IS strength, we note, in accord with Fig. 45, that the transitions appear 
more collective in the Ni and Sn drip-line isotopes than in that of Ca. The component 
with the largest amplitude {X^^ - Y^^ = 0.81) in the state at E ^ 0.713 MeV of ^^Ca is 
i^3si/2{£ = 0.355 MeV, vl = 0.65) ® p2d5/2{S = 0.618 MeV, = 0.78), and that with 
the second largest (X^^ — Y^^^ = 0.16) is {z/2(i5/2}^. This is still a predominantly single 
two-quasiparticle excitation. By contrast, the components of the state with at E = 0.096 
MeV in ^^Ni with the largest amplitudes are 

{ulhn/2{e = 1.415 MeY,vl = 0.27)}', - Y^ = 0.35, 
{vlg,l2{E = 1.350 MeV, vl = 0.59)}^ - Yl^ = 0.30, 

7r2p3/2(^: = 2.352 MeV, = 0) (g) i^lhi^i^ = 2.222 MeV, vl = 1), - Yl, = 0.12. 

While still not very collective, this state is a little more so than its counterpart in Ca. (^^Ni 
may be in the spherical-deformed transition region.) 

In ^''^Sn, the component of the state at E = 3.228 MeV with the largest amplitude is 
(Fig. 45) u2gg/2{S = 1.477 MeV, vl = 0)^iylin/2{S = 2.475 MeV, vl = 1.00), (X^^ -F^l = 



Ni has a small IV peak at the IS peak. 



36 



0.79) and that with the next largest is 7r2d5/2{S = 3.654 MeV, vf^ = 0) TJ^l99/2{S = 2.276 
MeV, vl = 1), {Xl^ - Y^^ = 0.09). This excitation is mainly a single neutron particle and 
neutron hole. Still, the proton contribution is much larger than that in other channels at 
the neutron drip line. This excited state has 23.7 % of the total IS strength and 8.6 % of 
the EWSR. 

We see both in this channel and others that the large transition strength near the neutron 
drip line does not always indicate the coherent contributions of many two-quasiparticle 
excitations (i.e., collectivity). Yet the transition densities make these strong low- lying states, 
particularly in Ni and Sn, look very much like the surface vibrations that characterize 
nuclei closer to stability. To get a better handle on the degree to which the states we 
examine here can be called vibrations, we show in Fig. 47 the energy, transition probability 
B{E2; 0"^ — > 2+), and largest — for neutrons in the lowest-energy 2+ states of Ca, 
Ni, and Sn. The Ca panel shows, first of all, that we do not reproduce experimental data 
very well except in a few nuclei, and second, that the states are dominated by a single 
two-quasiparticle excitation for < 34. For 34 < < 40 something peculiar happens; 
when the energy increases, so docs the B{E2;0'^ 2"*"), and vice versa, with the largest 
— Y^^, a measure of collectivity, varying irregularly. In nuclei with collective vibrations, 
transition strength is usually correlated inversely with energy [37, 38]; we see the opposite 
here. Near the neutron drip line of Ni, the classic signs of increasing quadrupole collectivity 
— decreasing energy, rising B{E2; 0+ 2+), and decreasing maximum X'^^ — Y"^^ — are all 
present. The Sn isotopes look much more familiar, exhibiting the usual trends, except at the 
magic numbers. And the states close to the drip line with 96 < A?" < 124 are as classically 
collective as (or more so than) those in more stable isotopes with 56 < < 68. 

Finally we show the interaction dependence of the collectivity measures for the lowest 
2"*" states in Sn (the easiest chain to understand) in Fig. 48. In the mid-shell region, the 
energies are low, the B{E2; 0"*" 2+)'s are high, and the largest amplitudes low, as expected 
from systematics in stable nuclei in regions of large dynamical deformation. The values of 
these quantities depend significantly on the interaction, however. Ca and Ni show similar 
behavior but with more irregularities. 
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FIG. 45: Transition densities to the excited states in the low-energy peak of the IS 2^ strength 
functions in the neutron-drip-line nuclei (wth SkM*) ^^Ca, ^^Ni, and ^^^Sn. Both proton and 
neutron curves are reduced by a factor of 5 in ^^Ni to fit on the figure.. 
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FIG. 46: Transition densities to the excited states in the low-energy peak of the IV 2+ strength 
functions (SkM*) in ^^^Sn and ''^Ni. 
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FIG. 47: iV-dependence of the energy (upper), B{E2;0^ —>■ 2+) (middle), and largest X'^^, - Y^^, 
(lower) of the lowest 2+ states in Ca, Ni, and Sn, with SkM*. Experimental data are also shown 
where available. 



39 




40 



V. CONCLUSION 



The feature that leaps out of our calculations is the enhancement of low-energy strength 
as the neutron drip line is approached. It is present in all multipoles (and is most pronounced 
in the IS 1~ channel, where the transition operator contains a factor r^), with all Skyrme 
interactions, and in every isotopic chain, beginning just after a closed neutron shell. Inter- 
estingly, the recent measurement of an IV 1~ strength function at GSI shows a low-energy 
peak before a closed shell, in ^^°Sn. This may be due to correlations not included in the 
QRPA. 

The transition densities to the low-lying states reveal some interesting features. In the 
O"*" and 1~ channels, the low-energy peaks are usually created essentially entirely by a single 
two-neutron-quasiparticle configuration. In the 0+ channel, where the total number of two- 
quasiparticle configurations is restricted, even small admixtures of two-proton-quasiparticle 
states are absent; the proton transition density is nearly zero. Protons play a noticeably 
larger role in the transition density of the 1~ channel despite the lack of any two-proton- 
quasiparticle configurations with sizable amplitude. In both these channels, the transition 
densities have a peak around the rms radius and very long tail outside the radius, which 
allows single quasiparticle configuration to have such large strength. 

The transition densities to 2"*" states arc different; they have no real nodes and a proton 
component that is of the same order as the neutron component. The transition densities 
near the neutron-drip line resemble those of surface vibrations, and in many cases the cor- 
responding states clearly are vibrational. But the usual systematics of collective vibrations 
is not always present; energies and E2 strengths are sometimes correlated in the wrong way. 

Some of the most significant differences in the predictions of the three Skyrme interac- 
tions we used are in static properties: the locations of the neutron drip line and quadrupole- 
deformed regions. But some dynamical properties also depend significantly on the interac- 
tion, most notably (and easiest to measure) the energies of and transition-strengths to the 
lowest 2+ states near regions of deformation and the total energy- weighted IV strength. 

This paper, in a way, is a catalog of predicted dynamic properties of nuclei near the 
drip line. For the most part we have not offered simple explanations for the anomalies that 
appear. That, we hope, will be the subject of a future paper. 
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APPENDIX: TRANSITION DENSITY 

The transition density is defined by 

NorZ 

1=1 
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where |0) and \k) denote the ground and excited states of the system, and in the right- 
hand side is an operator acting on the nucleon i. The transition density is evaluated in the 
QRPA (we omit q) as 

K<K' 

-Ykk' [^^K,{r)ipK{r)vK'UK - ip]iipK'{r)vKUK'') } , (A.2) 

where X'^^, and Y^j^, are the forward and backward amphtudes of the QRPA sohition k 
(see Appendix A of Ref. [28]), and ipxir) stands for a single-particle canonical- basis wave 
function. The uv factors associated with the basis are represented by uk and vk- Assuming 
spherical symmetry for the system and using 

^|;K{r) = i?^(r)F,,,(f))(//,^a|jm)|a), (A.3) 

ma 

with K = (fim) = {nljm) a set of quantum numbers, Ruir) a radial wave function, and \(t) 
a spin wave function, one can derive 



Ptr(r;/c) = ^ V(2i + l)(2/ + l)i^^Kr)i?» 

( -1/2 1/2 ) + 

where jj! — {n'j'V), and denotes the angular momentum of the state k. X^^, and Yj^^, are 
the forward and backward amplitudes in the spherical representation. 
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